Cell-attribute aware community detection improves differential abundance testing from single-cell RNA-Seq data.

Variations of cell-type proportions within tissues could be informative of biological aging and disease risk. Single-cell RNA-sequencing offers the opportunity to detect such differential abundance patterns, yet this task can be statistically challenging due to the noise in single-cell data, inter-sample variability and because such patterns are often of small effect size. Here we present a differential abundance testing paradigm called ELVAR that uses cell attribute aware clustering when inferring differentially enriched communities within the single-cell manifold. Using simulated and real single-cell and single-nucleus RNA-Seq datasets, we benchmark ELVAR against an analogous algorithm that uses Louvain for clustering, as well as local neighborhood-based methods, demonstrating that ELVAR improves the sensitivity to detect cell-type composition shifts in relation to aging, precancerous states and Covid-19 phenotypes. In effect, leveraging cell attribute information when inferring cell communities can denoise single-cell data, avoid the need for batch correction and help retrieve more robust cell states for subsequent differential abundance testing. ELVAR is available as an open-source R-package.

Bioinformatics and system biology approach to identify the influences among COVID-19, ARDS and sepsis.

Background Severe coronavirus disease 2019 (COVID -19) has led to severe pneumonia or acute respiratory distress syndrome (ARDS) worldwide. we have noted that many critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. The molecular mechanisms that underlie COVID-19, ARDS and sepsis are not well understood. The objectives of this study were to analyze potential molecular mechanisms and identify potential drugs for the treatment of COVID-19, ARDS and sepsis using bioinformatics and a systems biology approach. Methods Three RNA-seq datasets (GSE171110, GSE76293 and GSE137342) from Gene Expression Omnibus (GEO) were employed to detect mutual differentially expressed genes (DEGs) for the patients with the COVID-19, ARDS and sepsis for functional enrichment, pathway analysis, and candidate drugs analysis. Results We obtained 110 common DEGs among COVID-19, ARDS and sepsis. ARG1, FCGR1A, MPO, and TLR5 are the most influential hub genes. The infection and immune-related pathways and functions are the main pathways and molecular functions of these three diseases. FOXC1, YY1, GATA2, FOXL, STAT1 and STAT3 are important TFs for COVID-19. mir-335-5p, miR-335-5p and hsa-mir-26a-5p were associated with COVID-19. Finally, the hub genes retrieved from the DSigDB database indicate multiple drug molecules and drug-targets interaction. Conclusion We performed a functional analysis under ontology terms and pathway analysis and found some common associations among COVID-19, ARDS and sepsis. Transcription factors-genes interaction, protein-drug interactions, and DEGs-miRNAs coregulatory network with common DEGs were also identified on the datasets. We believe that the candidate drugs obtained in this study may contribute to the effective treatment of COVID-19.

Using bioinformatics and systems biology methods to identify the mechanism of interaction between COVID-19 and nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is considered a risk factor for severe COVID-19, but the mechanism remains unknown. This study used bioinformatics to help define the relationship between these diseases. The GSE147507 (COVID-19), GSE126848 (NAFLD), and GSE63067 (NAFLD-2) datasets were screened using the Gene Expression Omnibus. Common differentially expressed genes were then identified using a Venn diagram. Gene ontology analysis and KEGG pathway enrichment were performed on the differentially expressed genes. A protein-protein interaction network was also constructed using the STRING platform, and key genes were identified using the Cytoscape plugin. GES63067 was selected for validation of the results. Analysis of ferroptosis gene expression during the development of the 2 diseases and prediction of their upstream miRNAs and lncRNAs. In addition, transcription factors (TFs) and miRNAs related to key genes were identified. Effective drugs that act on target genes were found in the DSigDB. The GSE147507 and GSE126848 datasets were crossed to obtain 28 co-regulated genes, 22 gene ontology terms, 3 KEGG pathways, and 10 key genes. NAFLD may affect COVID-19 progression through immune function and inflammatory signaling pathways. CYBB was predicted to be a differential ferroptosis gene associated with 2 diseases, and the CYBB-hsa-miR-196a/b-5p-TUG1 regulatory axis was identified. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Eckol, sulfinpyrazone, and phenylbutazone) were considered as target drugs for Patients with COVID-19 and NAFLD. This study identified key gene and defined molecular mechanisms associated with the progression of COVID-19 and NAFLD. COVID-19 and NAFLD progression may regulate ferroptosis through the CYBB-hsa-miR-196a/b-5p-TUG1 axis. This study provides additional drug options for the treatment of COVID-19 combined with NAFLD disease.

Identifying key genes related to inflammasome in severe COVID-19 patients based on a joint model with random forest and artificial neural network.

Background: The coronavirus disease 2019 (COVID-19) has been spreading astonishingly and caused catastrophic losses worldwide. The high mortality of severe COVID-19 patients is an serious problem that needs to be solved urgently. However, the biomarkers and fundamental pathological mechanisms of severe COVID-19 are poorly understood. The aims of this study was to explore key genes related to inflammasome in severe COVID-19 and their potential molecular mechanisms using random forest and artificial neural network modeling. Methods: Differentially expressed genes (DEGs) in severe COVID-19 were screened from GSE151764 and GSE183533 via comprehensive transcriptome Meta-analysis. Protein-protein interaction (PPI) networks and functional analyses were conducted to identify molecular mechanisms related to DEGs or DEGs associated with inflammasome (IADEGs), respectively. Five the most important IADEGs in severe COVID-19 were explored using random forest. Then, we put these five IADEGs into an artificial neural network to construct a novel diagnostic model for severe COVID-19 and verified its diagnostic efficacy in GSE205099. Results: Using combining P value < 0.05, we obtained 192 DEGs, 40 of which are IADEGs. The GO enrichment analysis results indicated that 192 DEGs were mainly involved in T cell activation, MHC protein complex and immune receptor activity. The KEGG enrichment analysis results indicated that 192 GEGs were mainly involved in Th17 cell differentiation, IL-17 signaling pathway, mTOR signaling pathway and NOD-like receptor signaling pathway. In addition, the top GO terms of 40 IADEGs were involved in T cell activation, immune response-activating signal transduction, external side of plasma membrane and phosphatase binding. The KEGG enrichment analysis results indicated that IADEGs were mainly involved in FoxO signaling pathway, Toll-like receptor, JAK-STAT signaling pathway and Apoptosis. Then, five important IADEGs (AXL, MKI67, CDKN3, BCL2 and PTGS2) for severe COVID-19 were screened by random forest analysis. By building an artificial neural network model, we found that the AUC values of 5 important IADEGs were 0.972 and 0.844 in the train group (GSE151764 and GSE183533) and test group (GSE205099), respectively. Conclusion: The five genes related to inflammasome, including AXL, MKI67, CDKN3, BCL2 and PTGS2, are important for severe COVID-19 patients, and these molecules are related to the activation of NLRP3 inflammasome. Furthermore, AXL, MKI67, CDKN3, BCL2 and PTGS2 as a marker combination could be used as potential markers to identify severe COVID-19 patients.

Reconstruction of the cell pseudo-space from single-cell RNA sequencing data with scSpace.

Tissues are highly complicated with spatial heterogeneity in gene expression. However, the cutting-edge single-cell RNA-seq technology eliminates the spatial information of individual cells, which contributes to the characterization of cell identities. Herein, we propose single-cell spatial position associated co-embeddings (scSpace), an integrative method to identify spatially variable cell subpopulations by reconstructing cells onto a pseudo-space with spatial transcriptome references (Visium, STARmap, Slide-seq, etc.). We benchmark scSpace with both simulated and biological datasets, and demonstrate that scSpace can accurately and robustly identify spatially variated cell subpopulations. When employed to reconstruct the spatial architectures of complex tissue such as the brain cortex, the small intestinal villus, the liver lobule, the kidney, the embryonic heart, and others, scSpace shows promising performance on revealing the pairwise cellular spatial association within single-cell data. The application of scSpace in melanoma and COVID-19 exhibits a broad prospect in the discovery of spatial therapeutic markers.

Identification of biomarkers, pathways, and potential therapeutic targets for heart failure using next-generation sequencing data and bioinformatics analysis.

BACKGROUND: Heart failure (HF) is the most common cardiovascular diseases and the leading cause of cardiovascular diseases related deaths. Increasing molecular targets have been discovered for HF prognosis and therapy. However, there is still an urgent need to identify novel biomarkers. Therefore, we evaluated biomarkers that might aid the diagnosis and treatment of HF. METHODS: We searched next-generation sequencing (NGS) dataset (GSE161472) and identified differentially expressed genes (DEGs) by comparing 47 HF samples and 37 normal control samples using limma in R package. Gene ontology (GO) and pathway enrichment analyses of the DEGs were performed using the g: Profiler database. The protein-protein interaction (PPI) network was plotted with Human Integrated Protein-Protein Interaction rEference (HiPPIE) and visualized using Cytoscape. Module analysis of the PPI network was done using PEWCC1. Then, miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed by Cytoscape software. Finally, we performed receiver operating characteristic (ROC) curve analysis to predict the diagnostic effectiveness of the hub genes. RESULTS: A total of 930 DEGs, 464 upregulated genes and 466 downregulated genes, were identified in HF. GO and REACTOME pathway enrichment results showed that DEGs mainly enriched in localization, small molecule metabolic process, SARS-CoV infections, and the citric acid tricarboxylic acid (TCA) cycle and respiratory electron transport. After combining the results of the PPI network miRNA-hub gene regulatory network and TF-hub gene regulatory network, 10 hub genes were selected, including heat shock protein 90 alpha family class A member 1 (HSP90AA1), arrestin beta 2 (ARRB2), myosin heavy chain 9 (MYH9), heat shock protein 90 alpha family class B member 1 (HSP90AB1), filamin A (FLNA), epidermal growth factor receptor (EGFR), phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1), cullin 4A (CUL4A), YEATS domain containing 4 (YEATS4), and lysine acetyltransferase 2B (KAT2B). CONCLUSIONS: This discovery-driven study might be useful to provide a novel insight into the diagnosis and treatment of HF. However, more experiments are needed in the future to investigate the functional roles of these genes in HF.

A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology.

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

Gene function and cell surface protein association analysis based on single-cell multiomics data.

Single-cell transcriptomics provides researchers with a powerful tool to resolve the transcriptome heterogeneity of individual cells. However, this method falls short in revealing cellular heterogeneity at the protein level. Previous single-cell multiomics studies have focused on data integration rather than exploiting the full potential of multiomics data. Here we introduce a new analysis framework, gene function and protein association (GFPA), that mines reliable associations between gene function and cell surface protein from single-cell multimodal data. Applying GFPA to human peripheral blood mononuclear cells (PBMCs), we observe an association of epithelial mesenchymal transition (EMT) with the CD99 protein in CD4 T cells, which is consistent with previous findings. Our results show that GFPA is reliable across multiple cell subtypes and PBMC samples. The GFPA python packages and detailed tutorials are freely available at https://github.com/studentiz/GFPA.

Identification of key differentially expressed genes in SARS-CoV-2 using RNA-seq analysis with a systems biology approach.

COVID-19 is associated with dysregulation of several genes and signaling pathways. Based on the importance of expression profiling in identification of the pathogenesis of COVID-19 and proposing novel therapies for this disorder, we have employed an in silico approach to find differentially expressed genes between COVID-19 patients and healthy controls and their relevance with cellular functions and signaling pathways. We obtained 630 DEmRNAs, including 486 down-regulated DEGs (such as CCL3 and RSAD2) and 144 up-regulated DEGs (such as RHO and IQCA1L), and 15 DElncRNAs, including 9 down-regulated DElncRNAs (such as PELATON and LINC01506) and 6 up-regulated DElncRNAs (such as AJUBA-DT and FALEC). The PPI network of DEGs showed the presence of a number immune-related genes such as those coding for HLA molecules and interferon regulatory factors. Taken together, these results highlight the importance of immune-related genes and pathways in the pathogenesis of COVID-19 and suggest novel targets for treatment of this disorder.

The molecular mechanism of cardiac injury in SARS-CoV-2 infection: Focus on mitochondrial dysfunction.

BACKGROUND: Coronavirus disease 2019(COVID-19) caused a large number of infections worldwide. Although some patients recovered from the disease, some of the other problems that accompanied it, such as cardiac injury, could affect the patient's subsequent quality of life and prognosis. OBJECTIVES: To clarify the molecular mechanism of cardiac injury in SARS-CoV-2 Infection. METHODS: The RNA-Seq dataset (GSE184715) comparing expression profiling of Mock human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) and SARS-CoV-2-infected hiPSC-CMs was downloaded from Gene Expression Omnibus (GEO). Differentially expressed genes(DEGs) were performed by the R software. Degs were analyzed by enrichment analysis to clarify the affected pathways. Hub genes were screened out by a PPI network constructed from Degs. Finally, Connectivity Map was used to screen for the treatment of COVID-19 induced cardiac injury. RESULTS: 2705 differentially expressed genes were identified. Enrichment analysis confirmed that mitochondrial dysfunction was caused by SARS-CoV-2, meanwhile, cardiac muscle contraction was suppressed and NF-κB was activated. Based on the PPI network, 15 hub genes were identified. These 15 down-regulated hub genes were mainly involved in the reduced activity of complexes in the mitochondrial respiratory chain associated with mitochondrial dysfunction. Moreover, 5 candidate drugs were identified to treat cardiac injury. CONCLUSION: In conclusion, SARS-CoV-2 infection of cardiomyocytes causes mitochondrial dysfunction, including reduced mitochondrial respiratory chain complex activity and decreased ATP synthesis, leading to cardiomyocyte apoptosis, while the activated NF-κB also induced cytokine storms, ultimately resulting in cardiac injury.

A comprehensive analysis of gene expression profiling data in COVID-19 patients for discovery of specific and differential blood biomarker signatures.

COVID-19 is a newly recognized illness with a predominantly respiratory presentation. Although initial analyses have identified groups of candidate gene biomarkers for the diagnosis of COVID-19, they have yet to identify clinically applicable biomarkers, so we need disease-specific diagnostic biomarkers in biofluid and differential diagnosis in comparison with other infectious diseases. This can further increase knowledge of pathogenesis and help guide treatment. Eight transcriptomic profiles of COVID-19 infected versus control samples from peripheral blood (PB), lung tissue, nasopharyngeal swab and bronchoalveolar lavage fluid (BALF) were considered. In order to find COVID-19 potential Specific Blood Differentially expressed genes (SpeBDs), we implemented a strategy based on finding shared pathways of peripheral blood and the most involved tissues in COVID-19 patients. This step was performed to filter blood DEGs with a role in the shared pathways. Furthermore, nine datasets of the three types of Influenza (H1N1, H3N2, and B) were used for the second step. Potential Differential Blood DEGs of COVID-19 versus Influenza (DifBDs) were found by extracting DEGs involved in only enriched pathways by SpeBDs and not by Influenza DEGs. Then in the third step, a machine learning method (a wrapper feature selection approach supervised by four classifiers of k-NN, Random Forest, SVM, Naïve Bayes) was utilized to narrow down the number of SpeBDs and DifBDs and find the most predictive combination of them to select COVID-19 potential Specific Blood Biomarker Signatures (SpeBBSs) and COVID-19 versus influenza Differential Blood Biomarker Signatures (DifBBSs), respectively. After that, models based on SpeBBSs and DifBBSs and the corresponding algorithms were built to assess their performance on an external dataset. Among all the extracted DEGs from the PB dataset (from common PB pathways with BALF, Lung and Swab), 108 unique SpeBD were obtained. Feature selection using Random Forest outperformed its counterparts and selected IGKC, IGLV3-16 and SRP9 among SpeBDs as SpeBBSs. Validation of the constructed model based on these genes and Random Forest on an external dataset resulted in 93.09% Accuracy. Eighty-three pathways enriched by SpeBDs and not by any of the influenza strains were identified, including 87 DifBDs. Using feature selection by Naive Bayes classifier on DifBDs, FMNL2, IGHV3-23, IGLV2-11 and RPL31 were selected as the most predictable DifBBSs. The constructed model based on these genes and Naive Bayes on an external dataset was validated with 87.2% accuracy. Our study identified several candidate blood biomarkers for a potential specific and differential diagnosis of COVID-19. The proposed biomarkers could be valuable targets for practical investigations to validate their potential.

Application of Deep Learning on Single-cell RNA Sequencing Data Analysis: A Review.

Single-cell RNA sequencing (scRNA-seq) has become a routinely used technique to quantify the gene expression profile of thousands of single cells simultaneously. Analysis of scRNA-seq data plays an important role in the study of cell states and phenotypes, and has helped elucidate biological processes, such as those occurring during the development of complex organisms, and improved our understanding of disease states, such as cancer, diabetes, and coronavirus disease 2019 (COVID-19). Deep learning, a recent advance of artificial intelligence that has been used to address many problems involving large datasets, has also emerged as a promising tool for scRNA-seq data analysis, as it has a capacity to extract informative and compact features from noisy, heterogeneous, and high-dimensional scRNA-seq data to improve downstream analysis. The present review aims at surveying recently developed deep learning techniques in scRNA-seq data analysis, identifying key steps within the scRNA-seq data analysis pipeline that have been advanced by deep learning, and explaining the benefits of deep learning over more conventional analytic tools. Finally, we summarize the challenges in current deep learning approaches faced within scRNA-seq data and discuss potential directions for improvements in deep learning algorithms for scRNA-seq data analysis.

Identification of a key glioblastoma candidate gene, FUBP3, based on weighted gene co-expression network analysis.

BACKGROUND: Glioblastoma multiforme (GBM) is the most common aggressive malignant brain tumor. However, the molecular mechanism of glioblastoma formation is still poorly understood. To identify candidate genes that may be connected to glioma growth and development, weighted gene co-expression network analysis (WGCNA) was performed to construct a gene co-expression network between gene sets and clinical characteristics. We also explored the function of the key candidate gene. METHODS: Two GBM datasets were selected from GEO Datasets. The R language was used to identify differentially expressed genes. WGCNA was performed to construct a gene co-expression network in the GEO glioblastoma samples. A custom Venn diagram website was used to find the intersecting genes. The GEPIA website was applied for survival analysis to determine the significant gene, FUBP3. OS, DSS, and PFI analyses, based on the UCSC Cancer Genomics Browser, were performed to verify the significance of FUBP3. Immunohistochemistry was performed to evaluate the expression of FUBP3 in glioblastoma and adjacent normal tissue. KEGG and GO enrichment analyses were used to reveal possible functions of FUBP3. Microenvironment analysis was used to explore the relationship between FUBP3 and immune infiltration. Immunohistochemistry was performed to verify the results of the microenvironment analysis. RESULTS: GSE70231 and GSE108474 were selected from GEO Datasets, then 715 and 694 differentially expressed genes (DEGs) from GSE70231 and GSE108474, respectively, were identified. We then performed weighted gene co-expression network analysis (WGCNA) and identified the most downregulated gene modules of GSE70231 and GSE108474, and 659 and 3915 module genes from GSE70231 and GSE108474, respectively, were selected. Five intersection genes (FUBP3, DAD1, CLIC1, ABR, and DNM1) were calculated by Venn diagram. FUBP3 was then identified as the only significant gene by survival analysis using the GEPIA website. OS, DSS, and PFI analyses verified the significance of FUBP3. Immunohistochemical analysis revealed FUBP3 expression in GBM and adjacent normal tissue. KEGG and GO analyses uncovered the possible function of FUBP3 in GBM. Tumor microenvironment analysis showed that FUBP3 may be connected to immune infiltration, and immunohistochemistry identified a positive correlation between immune cells (CD4 + T cells, CD8 + T cells, and macrophages) and FUBP3. CONCLUSION: FUBP3 is associated with immune surveillance in GBM, indicating that it has a great impact on GBM development and progression. Therefore, interventions involving FUBP3 and its regulatory pathway may be a new approach for GBM treatment.

Transcriptomics for child and adolescent tuberculosis.

Tuberculosis (TB) in humans is caused by Mycobacterium tuberculosis (Mtb). It is estimated that 70 million children (<15 years) are currently infected with Mtb, with 1.2 million each year progressing to disease. Of these, a quarter die. The risk of progression from Mtb infection to disease and from disease to death is dependent on multiple pathogen and host factors. Age is a central component in all these transitions. The natural history of TB in children and adolescents is different to adults, leading to unique challenges in the development of diagnostics, therapeutics, and vaccines. The quantification of RNA transcripts in specific cells or in the peripheral blood, using high-throughput methods, such as microarray analysis or RNA-Sequencing, can shed light into the host immune response to Mtb during infection and disease, as well as understanding treatment response, disease severity, and vaccination, in a global hypothesis-free manner. Additionally, gene expression profiling can be used for biomarker discovery, to diagnose disease, predict future disease progression and to monitor response to treatment. Here, we review the role of transcriptomics in children and adolescents, focused mainly on work done in blood, to understand disease biology, and to discriminate disease states to assist clinical decision-making. In recent years, studies with a specific pediatric and adolescent focus have identified blood gene expression markers with diagnostic or prognostic potential that meet or exceed the current sensitivity and specificity targets for diagnostic tools. Diagnostic and prognostic gene expression signatures identified through high-throughput methods are currently being translated into diagnostic tests.

A comparative study of COVID-19 transcriptional signatures between clinical samples and preclinical cell models in the search for disease master regulators and drug repositioning candidates.

Coronavirus disease 2019 (COVID-19) is an acute viral disease with millions of cases worldwide. Although the number of daily new cases and deaths has been dropping, there is still a need for therapeutic alternatives to deal with severe cases. A promising strategy to prospect new therapeutic candidates is to investigate the regulatory mechanisms involved in COVID-19 progression using integrated transcriptomics approaches. In this work, we aimed to identify COVID-19 Master Regulators (MRs) using a series of publicly available gene expression datasets of lung tissue from patients which developed the severe form of the disease. We were able to identify a set of six potential COVID-19 MRs related to its severe form, namely TAL1, TEAD4, EPAS1, ATOH8, ERG, and ARNTL2. In addition, using the Connectivity Map drug repositioning approach, we identified 52 different drugs which could be used to revert the disease signature, thus being candidates for the design of novel clinical treatments. Furthermore, we compared the identified signature and drugs with the ones obtained from the analysis of nasopharyngeal swab samples from infected patients and preclinical cell models. This comparison showed significant similarities between them, although also revealing some limitations on the overlap between clinical and preclinical data in COVID-19, highlighting the need for careful selection of the best model for each disease stage.

Robust single-cell matching and multimodal analysis using shared and distinct features.

The ability to align individual cellular information from multiple experimental sources is fundamental for a systems-level understanding of biological processes. However, currently available tools are mainly designed for single-cell transcriptomics matching and integration, and generally rely on a large number of shared features across datasets for cell matching. This approach underperforms when applied to single-cell proteomic datasets due to the limited number of parameters simultaneously accessed and lack of shared markers across these experiments. Here, we introduce a cell-matching algorithm, matching with partial overlap (MARIO) that accounts for both shared and distinct features, while consisting of vital filtering steps to avoid suboptimal matching. MARIO accurately matches and integrates data from different single-cell proteomic and multimodal methods, including spatial techniques and has cross-species capabilities. MARIO robustly matched tissue macrophages identified from COVID-19 lung autopsies via codetection by indexing imaging to macrophages recovered from COVID-19 bronchoalveolar lavage fluid by cellular indexing of transcriptomes and epitopes by sequencing, revealing unique immune responses within the lung microenvironment of patients with COVID.

Discovering common pathogenetic processes between COVID-19 and sepsis by bioinformatics and system biology approach.

Corona Virus Disease 2019 (COVID-19), an acute respiratory infectious disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has spread rapidly worldwide, resulting in a pandemic with a high mortality rate. In clinical practice, we have noted that many critically ill or critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. In addition, it has been demonstrated that severe COVID-19 has some pathological similarities with sepsis, such as cytokine storm, hypercoagulable state after blood balance is disrupted and neutrophil dysfunction. Considering the parallels between COVID-19 and non-SARS-CoV-2 induced sepsis (hereafter referred to as sepsis), the aim of this study was to analyze the underlying molecular mechanisms between these two diseases by bioinformatics and a systems biology approach, providing new insights into the pathogenesis of COVID-19 and the development of new treatments. Specifically, the gene expression profiles of COVID-19 and sepsis patients were obtained from the Gene Expression Omnibus (GEO) database and compared to extract common differentially expressed genes (DEGs). Subsequently, common DEGs were used to investigate the genetic links between COVID-19 and sepsis. Based on enrichment analysis of common DEGs, many pathways closely related to inflammatory response were observed, such as Cytokine-cytokine receptor interaction pathway and NF-kappa B signaling pathway. In addition, protein-protein interaction networks and gene regulatory networks of common DEGs were constructed, and the analysis results showed that ITGAM may be a potential key biomarker base on regulatory analysis. Furthermore, a disease diagnostic model and risk prediction nomogram for COVID-19 were constructed using machine learning methods. Finally, potential therapeutic agents, including progesterone and emetine, were screened through drug-protein interaction networks and molecular docking simulations. We hope to provide new strategies for future research and treatment related to COVID-19 by elucidating the pathogenesis and genetic mechanisms between COVID-19 and sepsis.

Blood transcriptome analysis revealed the crosstalk between COVID-19 and HIV.

Background: The severe coronavirus disease 2019 (COVID-19) is an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has resulted in the most devastating pandemic in modern history. Human immunodeficiency virus (HIV) destroys immune system cells and weakens the body's ability to resist daily infections and diseases. Furthermore, HIV-infected individuals had double COVID-19 mortality risk and experienced worse COVID-related outcomes. However, the existing research still lacks the understanding of the molecular mechanism underlying crosstalk between COVID-19 and HIV. The aim of our work was to illustrate blood transcriptome crosstalk between COVID-19 and HIV and to provide potential drugs that might be useful for the treatment of HIV-infected COVID-19 patients. Methods: COVID-19 datasets (GSE171110 and GSE152418) were downloaded from Gene Expression Omnibus (GEO) database, including 54 whole-blood samples and 33 peripheral blood mononuclear cells samples, respectively. HIV dataset (GSE37250) was also obtained from GEO database, containing 537 whole-blood samples. Next, the "Deseq2" package was used to identify differentially expressed genes (DEGs) between COVID-19 datasets (GSE171110 and GSE152418) and the "limma" package was utilized to identify DEGs between HIV dataset (GSE37250). By intersecting these two DEG sets, we generated common DEGs for further analysis, containing Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) functional enrichment analysis, protein-protein interaction (PPI) analysis, transcription factor (TF) candidate identification, microRNAs (miRNAs) candidate identification and drug candidate identification. Results: In this study, a total of 3213 DEGs were identified from the merged COVID-19 dataset (GSE171110 and GSE152418), and 1718 DEGs were obtained from GSE37250 dataset. Then, we identified 394 common DEGs from the intersection of the DEGs in COVID-19 and HIV datasets. GO and KEGG enrichment analysis indicated that common DEGs were mainly gathered in chromosome-related and cell cycle-related signal pathways. Top ten hub genes (CCNA2, CCNB1, CDC20, TOP2A, AURKB, PLK1, BUB1B, KIF11, DLGAP5, RRM2) were ranked according to their scores, which were screened out using degree algorithm on the basis of common DEGs. Moreover, top ten drug candidates (LUCANTHONE, Dasatinib, etoposide, Enterolactone, troglitazone, testosterone, estradiol, calcitriol, resveratrol, tetradioxin) ranked by their P values were screened out, which maybe be beneficial for the treatment of HIV-infected COVID-19 patients. Conclusion: In this study, we provide potential molecular targets, signaling pathways, small molecular compounds, and promising biomarkers that contribute to worse COVID-19 prognosis in patients with HIV, which might contribute to precise diagnosis and treatment for HIV-infected COVID-19 patients.

In Silico Identification and Characterization of circRNAs During Host-Pathogen Interactions.

Circular RNAs (circRNAs) are a class of non-coding RNAs that are formed via back-splicing. These circRNAs are predominantly studied for their roles as regulators of various biological processes. Notably, emerging evidence demonstrates that host circRNAs can be differentially expressed (DE) upon infection with pathogens (e.g., influenza and coronaviruses), suggesting a role for circRNAs in regulating host innate immune responses. However, investigations on the role of circRNAs during pathogenic infections are limited by the knowledge and skills required to carry out the necessary bioinformatic analysis to identify DE circRNAs from RNA sequencing (RNA-seq) data. Bioinformatics prediction and identification of circRNAs is crucial before any verification, and functional studies using costly and time-consuming wet-lab techniques. To solve this issue, a step-by-step protocol of in silico prediction and characterization of circRNAs using RNA-seq data is provided in this manuscript. The protocol can be divided into four steps: 1) Prediction and quantification of DE circRNAs via the CIRIquant pipeline; 2) Annotation via circBase and characterization of DE circRNAs; 3) CircRNA-miRNA interaction prediction through Circr pipeline; 4) functional enrichment analysis of circRNA parental genes using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). This pipeline will be useful in driving future in vitro and in vivo research to further unravel the role of circRNAs in host-pathogen interactions.

Interaction-based transcriptome analysis via differential network inference.

Gene-based transcriptome analysis, such as differential expression analysis, can identify the key factors causing disease production, cell differentiation and other biological processes. However, this is not enough because basic life activities are mainly driven by the interactions between genes. Although there have been already many differential network inference methods for identifying the differential gene interactions, currently, most studies still only use the information of nodes in the network for downstream analyses. To investigate the insight into differential gene interactions, we should perform interaction-based transcriptome analysis (IBTA) instead of gene-based analysis after obtaining the differential networks. In this paper, we illustrated a workflow of IBTA by developing a Co-hub Differential Network inference (CDN) algorithm, and a novel interaction-based metric, pivot APC2. We confirmed the superior performance of CDN through simulation experiments compared with other popular differential network inference algorithms. Furthermore, three case studies are given using colorectal cancer, COVID-19 and triple-negative breast cancer datasets to demonstrate the ability of our interaction-based analytical process to uncover causative mechanisms.